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Abstract 

We calculate the 0(a s ) and 0{awfn1 /^w) corrections to the production of a single 
top quark via the weak process qq — > tb at the Fermilab Tevatron and the CERN Large 
Hadron Collider. An accurate calculation of the cross section is necessary in order to 
extract \Vtb\ from experiment. 



Figure 1: Single-top-quark production via qq — > i&. 

1 Introduction 

The recent discovery of the top quark [l[ has focused attention on top-quark physics. 
With the advent of accelerators able to produce copious numbers of top quarks, a comparison 
of the top quark's observed properties with those predicted by the Standard Model promises 
to be an important test of the model and may well provide insight into exciting new physics. 

In this paper we calculate the next-to-leading-order cross section for the weak process 
qq — > tb, which produces a single top quark via a virtual s-channel W boson (Fig. P [§, |J. 
The most important corrections to the 0{a 2 v ) leading-order cross section are the QCD 
correction of 0(a s ) and the Yukawa correction of 0(awm 2 /M^). The Yukawa correction, 
which arises from loops of Higgs bosons and the scalar components of virtual vector bosons, 
dominates the ordinary 0(aw) electroweak correction in the large m t limit. For the known 
value of the top-quark mass, m t = 175 ± 6 GeV, the Yukawa correction is expected to be at 
least as large as the ordinary electroweak correction. 

A precise theoretical calculation of the cross section for qq — > tb is necessary for a number 
of reasons. The cross section obviously determines the yield of single top quarks produced 
via this process. More importantly, the coupling of the top quark to the W boson in qq — > tb 
is proportional to the Cabibbo-Kobayashi-Maskawa (CKM) matrix element Vtb, one of the 
few Standard Model parameters not yet measured experimentally. If there are only three 
generations, unitarity of the CKM matrix implies that \Vtb\ must be very close to unity 
(.9988 < | Vtb | < .9995) ||. However, if there is a fourth generation, |Vtb| could be anything 
between (almost) zero and unity, depending on the amount of mixing between the third and 
fourth generations. Measurement of the qq —> tb cross section, coupled with an accurate 
theoretical calculation, may provide the best direct measurement of |Vtb| [§]■ Finally, in 
addition to being interesting in its own right, qq — > tb is a significant background to other 
processes, such as qq — > WH with H — > bb, where H is the Higgs boson ||. 

In some ways, qq — > tb is similar to the more-studied W-gluon fusion process (Fig. 0) ||. 
However, where that process involves a space-like W boson with q 2 < 0, the process qq — > tb 
proceeds via a time-like W boson with q 2 > (m t + mb) 2 . Thus these two processes, together 
with the decay of the top quark, t — > Wb (where the W boson has q 2 ~ Mjy), probe comple- 
mentary aspects of the top quark's weak charged current. The kinematic distributions of the 
final-state particles in the two processes also differ significantly. There is an additional jet 
present in W-gluon fusion, and the b quark is usually produced at low transverse momentum, 
while in qq — > tb, the b quark recoils against the t quark with high transverse momentum. 

At the Fermilab Tevatron (^/S = 2 TeV pp collider), the sum of the cross sections for 
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Figure 2: Single-top-quark production via W^-gluon fusion. 

qq — > tb and qq — > ib is roughly a factor of seven smaller than the dominant ti production 
cross section [^], and about a factor of two smaller than the PU-gluon- fusion cross section 
@. Nevertheless, a recent study indicates that with double b tagging, a signal is observable 
at the Tevatron with 2-3 fb _1 of integrated luminosity Unfortunately, even though the 
qq — > tb, ib cross section is larger at the CERN Large Hadron Collider (LHC, v^S* = 14 TeV 
pp collider), the signal will likely be obscured by backgrounds from the even larger ti and 
PU-gluon fusion processes, which are initiated by gluons ||. 

An important feature of qq — > tb is the accuracy with which the cross section can be 
calculated. The top-quark mass is much larger than Aqcd, so calculations are performed in a 
regime where perturbative QCD is very reliable. The correction to the initial state is identical 
to that occurring in the ordinary Drell-Yan process qq — > W* — > Eu (W* denotes a virtual W 
boson), which has been calculated to 0(a 2 s ) |J. Furthermore, by experimentally measuring 
qq — > W* —>■ £u, the initial quark- ant iquark flux can be constrained without recourse to 
perturbation theory^ This provides a check of the parton distribution functions, and allows 
the reduction of systematic errors. The parton distribution functions are not expected to be 
a large source of uncertainty, as the dominant contribution to the cross section comes from 
quark and antiquark distribution functions evaluated at relatively high values of x, where 
they are well known. There is little sensitivity to the less-well-known gluon distribution 
function, in contrast to the case of P^-gluon fusion. The final-state correction to the inclusive 
cross section is straightforward, and involves no collinear or infrared singularities. The QCD 
corrections to the initial and final states do not interfere at next-to-leading-order because 
the tb is in a color singlet if a gluon is emitted from the initial state, but a color octet if it is 
emitted from the final state. There is, however, interference at 0{a 2 s ) from the emission of 
two gluons. 

This paper is organized as follows. In Section 2 we present the 0(a s ) QCD corrections 
to both the initial and final states, and discuss their dependence on the renormalization 
and factorization scales. In Section 3 we present the 0{ay/m 2 / My/) Yukawa correction. 

1 Since the longitudinal momentum of the neutrino cannot be reconstructed, the q 2 of the W* cannot be 
determined, so qq —> W* — > Iv yields only a constraint on the quark- antiquark flux, rather than a direct 
measurement. 
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Figure 3: 0(a s ) correction to qq ~ * tb: (a)-(c) initial state, (d) final state. 



In Section 4 we present a summary of our results. We give an analytic expression for the 
Yukawa correction in an appendix. 



2 QCD correction 

The diagrams which contribute to the 0(a s ) correction to qq — > tb are shown in Fig. || 
As mentioned in the Introduction, the QCD corrections to the initial and final states do not 
interfere at 0(a s ). Therefore, we may consider the corrections to the initial and final states 
separately. To this end, we break up the process pp — > tb + X into the production of a virtual 
W boson of mass-squared q 2 , followed by its propagation and decay into tb. 

The production cross section of the virtual W boson is formally identical to that of the 
Drell-Yan process, to all orders in QCD. The modulus squared of the decay amplitude of the 
virtual W boson, integrated over the phase space of all final-state particles, is obtained by 
the application of Cutkosky's rules as twice the imaginary part of the self-energy of the 
W boson due to a tb loop, again to all orders in QCD. Furthermore, because the current to 
which the W boson couples in the initial state is conserved to all orders in QCD (for massless 
quarks), we need only consider the —g^ u term in the P^-boson propagator and self-energy. 
Thus we may write the differential cross section as 

dcr , ,7 ^\ / TT ,* , r s I m n(g 2 , m?, ml) /lX 

-fep ^a + X ) = a(0 - w + x) t( w_ (1) 

where II is the coefficient of the — g^ v term of the self-energy of a W boson with mass- 
squared q 2 . The total cross section is obtained by integrating over q 2 . This equation is valid 
to 0(a s ), but not beyond, because it neglects the interference between the QCD corrections 
to the initial and final states. 
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To demonstrate this procedure, we obtain the leading-order cross section for pp — > tb 
using 

a(pp^W*+X) = ^2 J dxi J dx 2 [qi (xx, /j lF )q j (x 2 , hf) + q~i(x!, [M F )qj(x 2 , /i f )] (2) 

x\Vv\ 2 ^5{ Xl z 2 S - f) 

where aw = g 2 /^ = V^G^M w /tt, S is the square of the total hadronic center-of-mass 
energy, q and q are the parton distribution functions, hf is the factorization scale, and the 
sum on % and j runs over all contributing quark-antiquark combinations. At leading order, 
the coefficient of the — g^ v term in the imaginary part of the l^-boson self-energy is 



Im Il(g 2 , m 2 , ml) 



a w \ l ' 2 \V tb \ 2 



ml + ml (ml — ml) 2 



2q 2 



2q 4 



(3) 



(4) 



2 (pp~^tb + X) = J2 / dx i / dx 2 [Qifa, fMF^fa, Pf) + qi(xi, fi F )qj(x 2 , Hf)] (5) 



where A is the triangle function associated with two-particle phase space, 

A = \(q 2 , m 2 , ml) = q 4 + m\ + m\ — 2q 2 m 2 t — 2q 2 ml — 2m\ml . 
Using Eq. (1), the differential cross section is thus 

da 
dq 



1, |2 f^AVW 
X| ^' 12(q 2 — M w ) 2 



1 - 



ml + ml 
2q 2 



m 



2\2 



2q 4 



5(xiX 2 S' - q 2 ) 



At leading order, the integration over q 2 to obtain the total cross section is trivial due to the 
delta function. At next-to-leading order, however, it is necessary to perform the integration 
numerically. 

The 0(a s ) corrections to the Drell-Yan process and the IV-boson self energy [JO, |T2 
were both calculated many years ago. We use the expression for a(pp — > W* + X) as given 
in Eqs. (9.5) and (12.3) of Ref. pj, and Im IT as derived from Eq. (3.3) of Ref. @fl We 
use m t = 175 GeV, m b = 5 GeV, M w = 80.33 GeV, \V tb \ — 1, — 1.16639 x 1(T 5 GeV" 2 , 
and a s as given by the parton distribution functions. 

The calculation of the initial-state correction includes divergences arising from collinear 
parton emission. These divergences cancel with corresponding divergences present in the 
QCD correction to the parton distribution functions. The finite terms remaining depend 
on the factorization scale /xp, both through the parton distribution functions and explicitly 
in the partonic cross section. The variation of the leading-order and next-to-leading-order 
cross sections with fip/y^q 2 , where yfq 2 is the mass of the virtual W boson,[] is shown 
in Fig. f| at both the Tevatron and the LHC. The leading-order cross section is calculated 

The exact correspondence between our notation and that of Ref. || is Im II = 37ra w |Vti,| 2 Im(Ilf + 11^). 
3 We have chosen to refer the scale \x F to the q 2 of the virtual W boson because this is the quantity which 
appears in the factorization logarithms. Thus the factorization scale hf varies when integrating over q 2 to 
obtain the total cross section. 
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Figure 4: Factorization-scale dependence of the leading-order (LO) and next-to-leading-order (NLO) cross 
sections for qq — > tb, tb at the Tevatron and the LHC. The NLO cross sections include only the initial-state 
QCD correction, and not the final-state correction. The LO cross sections are calculated with the CTEQ3L 
LO parton distribution functions, and the NLO cross sections with the CTEQ3M NLO parton distribution 
functions. 



with the CTEQ3L leading-order parton distribution functions, and the next-to-leading-order 



cross section with the CTEQ3M next-to-leading-order parton distribution functions [|14 
The leading-order cross section varies considerably with n?, while the next-to- leading-order 
cross section is appreciably less sensitive. The next-to-leading-order cross section shown 
in Fig. 4 contains only the initial-state correction. We see that for fi F = yfq 1 the initial- 
state correction is +36% at the Tevatron and +33% at the LHC|] In what follows, we set 

The cross section at next-to-leading order also depends on the renormalization scale, fiR, 
at which a s is evaluated. In Fig. [5] we show the next-to- leading-order cross section, including 
both initial- and final-state corrections, as a function of fin/y/q 1 , at both the Tevatron and 
the LHC. The dependence of the cross section on the renormalization scale first appears at 
next-to-leading order and is therefore mild. In what follows, we set fiR = \fq I . The final- 
state correction increases the cross section by +18% of the leading-order cross section at the 
Tevatron and +17% at the LHC. 

We show in Fig. ^ the leading-order and next-to-leading-order differential cross section as 
a function of the mass of the virtual W boson, ^fq 1 , at both the Tevatron and the LHC. Also 
shown are the separate 0(a s ) corrections from the initial and final states. These corrections 
have different shapes from the leading-order cross section, and from each other. In order 
to observe tb production experimentally, it is necessary to detect the b quark ||. Thus 

4 If both the leading-order and next-to-leading-order cross sections are calculated with the CTEQ3M 
next-to-leading-order parton distribution functions, the initial-state correction is +27% at the Tevatron and 
+15% at the LHC. Thus +9% of the initial-state correction at the Tevatron, and +18% at the LHC, is 
due to the increase in the leading-order cross section when it is calculated with next-to-leading-order parton 
distribution functions. 
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Figure 5: Renormalization-scale dependence of the leading-order (LO) and next-to-leading-order (NLO) 
cross sections for qq — > tb,ib at the Tevatron and the LHC. The NLO cross sections include both the 
initial-state and final-state correction. The LO cross sections are calculated with the CTEQ3L LO parton 
distribution functions, and the NLO cross sections with the CTEQ3M NLO parton distribution functions. 



the measured cross section will exclude some region near threshold, where the b quark does 
not have sufficient transverse momentum to be detected with high efficiency. Therefore the 
measured cross section, as well as the QCD correction, will depend on the acceptance for 
the b quark. 

If the top and bottom quarks were stable, they would form quarkonium bound states 
just below threshold |T5[. We estimate the distance below threshold that the ground state 



would occur, by analogy with the hydrogen atom, to be E « (4a s /3) 2 mb/2 « 50 MeV.Q The 
formation time of the ground state is approximately 1/E. This is much greater than the 
top-quark lifetime, IY ^ (1.5 GeV) -1 , so there is not sufficient time for quarkonium bound 



states to form 16 



Because the top-quark width is small compared to its mass, interference between the 
corrections to production and decay amplitudes has a negligible effect, of order a s T t /m t , on 
the total cross section |I7| |. This interference also has a negligible effect on differential cross 
sections, such as the invariant-mass distribution of the decay products of the top quark JL8 



Our final results for the cross section and uncertainty will be presented in Section 4. 



3 Yukawa correction 

The diagrams which contribute to the 0(awmt / m w) Yukawa correction to qq — > tb are 
shown in Fig. [|]. The dashed lines represent the Higgs boson and the unphysical scalar W 
and Z bosons associated with the Higgs field (in the gauge). The effect of a top-quark 

5 Here is the approximate reduced mass of the system, and Gp = 4/3 is the usual SU(3) group theory 
factor associated with the fundamental representation. 
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Figure 6: Differential cross section for qq — ► tb,ib versus the mass of the virtual s-channel W boson, at 
the Tevatron and the LHC. Both the leading-order (LO) and next-to-leading order (NLO) cross sections are 
shown, as well as the separate contributions from the initial-state (IS) and final-state (FS) corrections. The 
LO cross sections are calculated with the CTEQ3L LO parton distribution functions, and the NLO cross 
sections with the CTEQ3M NLO parton distribution functions. 





Figure 7: 0(aw r m'^/M^) corrections to qq — > tb. The dashed lines represent the Higgs boson and the 
unphysical scalar W and Z bosons (in gauge) : (a) wavefunction renormalization, (b) vertex correction. 



loop in the P^-boson propagator, which might be expected to contribute a term of Yukawa 
strength, is absorbed by the renormalized weak coupling constant, which we express in terms 
of Gfj,, the Fermi constant measured in muon decay (aw = fl ,2 /4vr = \/2G \iM w j V) . We use 
standard Feynman integral techniques with dimensional regularization to calculate the loop 
diagrams [TI|, and work in the approximation where the bottom quark is massless. Our other 



parameters are m t = 175 GeV, M w = 80.33 GeV, \V tb \ = 1, and G M = 1.16639 x 1(T 5 GeV 
In the m& = approximation, the matrix element of the tb current may be written as 



-2 



iu(p t )T^v(p b ) 



-igT 



2V2 



+ 




(6) 



I u(p t 
8V2tt 2 1 



(1 - 7>(Pb) 



where pt and p b are the outgoing four-momenta of the t and b quarks, respectively; the form 
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Figure 8: Fractional change in the total cross section for qq 
Higgs-boson mass at the Tevatron and the LHC. 



tb, tb due to the Yukawa correction vs. the 



factors F\ and F2 are functions of q 2 



SU(3) matrix [Tr(T A T B ^ 



(pt + Pb) 2 and the Higgs-boson mass; and T A is an 



¥ AB \ 



The fractional change in the differential cross section as 



a function of the q of the virtual W boson is 



Adcry I dy/q 2 



m 2 G^ 



da LO /dy/q 



8y/2- 



2\2 



2F 1 (q 2 ) + F 2 (q 2 



(q - m t) 



q 4 - 



1 9 2 



(7) 



Analytic expressions for the form factors F\ and F2 are given in an appendix. 

The fractional change in the total cross section, Aay/uLO, is plotted in Fig. || vs. the 
Higgs-boson mass, Mh, at both the Tevatron and the LHC. For values of Mh between 60 
GeV and 1 TeV, the absolute value of the Yukawa correction is never more than one percent 
of the leading-order cross section. Thus the Yukawa correction is negligible for this process, 
as has also been found to be the case for tt production [^D], |2T|. Since W-gluon fusion also 



involves the tb weak charged current, our calculation suggests that the Yukawa correction 
to that process is also negligible. As previously mentioned, the ordinary weak correction is 
expected to be comparable to the Yukawa correction, so it too should be negligible.^ The 
Yukawa correction could potentially be much larger in models with enhanced couplings of 
Higgs bosons to top or bottom quarks pO, E2 . 



4 Conclusions 

The cross section for qq — > tb, tb at both the Tevatron and the LHC is given in Table 
|TJ. The leading-order cross section, next-to-leading-order cross section including only the 
initial-state QCD correction, and the full next-to-leading order cross section are given. The 

6 The complete calculation of the ordinary weak correction would require a set of parton distribution 
functions which are extracted with weak corrections included. Such a set is not available at this time. 
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m t = 175 GeV, [ir = [If = Vy 1 


CTEQ3L,3M 


MRS (A') 


MRS(G) 


Tevatron 
VS = 2 TeV 




.578 


.601 


.602 


CNLO{IS) 


.789 


.766 


.758 


O'NLO 


.894 


.868 


.860 


cnlo (avg) = .88 ± .05 pb 


LHC 

= 14 TeV 




6.76 


7.83 


7.81 


VNLO(IS) 


9.02 


9.01 


9.03 


&NLO 


10.19 


10.17 


10.21 


(Tnlo (avg) = 10.2 ± 0.6 pb 



Table 1 : Leading-order (LO) and next-to- leading-order (NLO) cross sections (pb) for qq tb, tb at the 
Tevatron and the LHC for three different sets of parton distribution functions (PDFs). The NLO cross 
section including only the initial state (IS) correction is also given. The CTEQ LO cross section is computed 
with the CTEQ3L LO PDFs; all other cross sections are computed with NLO PDFs. The final NLO cross 
section is the average of the CTEQ3M and MRS(A') cross sections, with an uncertainty of ±6%, as discussed 
in the text. 



factorization and renormalization scales are both set equal to yip, the mass of the virtual 
W boson. We give results for three different sets of next-to-leading-order parton distribu- 
tion functions: CTEQ3M 0, MRS (A'), and MRS(G) H|] The QCD correction to the 
cross section is significant: about +54% at the Tevatron and +50% at the LHC, with the 
leading-order cross section evaluated with leading-order parton distribution functions, and 
the next-to-leading-order cross section evaluated with next-to-leading-order parton distribu- 
tion functions.^ The size of the 0(a s ) correction improves the outlook for observation of 
this process in Run II at the Tevatron. 

As shown in Fig. |], varying the factorization scale between one half and twice y/q 1 changes 
the cross section by only ±2%. Varying the renormalization scale over this same range yields 
a similar change in the cross section, as shown in Fig. [5]. Using these results to estimate 
the contribution from higher-order QCD corrections, we conclude that the uncertainty in 
the cross section is at the level of ±4%. This conclusion is supported by the known next-to- 
next-to-leading-order correction to the Drell-Yan process, which is about 2% (in the modified 
minimal subtraction (MS) scheme) ||. 

It is difficult to reliably ascertain the uncertainty in the cross section from the parton 
distribution functions at this time. The small difference in the next-to-leading-order cross 
sections using MRS (A ) and MRS(G) supports the contention that the calculation is in- 
sensitive to the gluon distribution function. The difference between the cross section using 
CTEQ3M and MRS (A ) suggests that the uncertainty in the cross section from the parton 
distribution functions is on the order of ±2%. However, since each set of parton distribution 
functions represents the best fit to some set of data, the uncertainty is certainly larger than 
this. Therefore, we assign an uncertainty of ±4% from the parton distribution functions. 

7 The leading-order CTEQ cross section is calculated with the CTEQ3L leading-order parton distribution 
functions. 

8 If next-to-leading-order parton distribution functions are used at both leading and next-to-leading order, 
the correction is about +45% at the Tevatron and +32% at the LHC. 



9 



Top-quark mass dependence 



Top-quark mass dependence 




m t (GeV) m t (GeV) 

Figure 9: Next- to- leading order cross section for qq — ► tb, tb as a function of the top-quark mass. 



For our final estimate, we average the next-to-leading-order cross sections using the 
CTEQ3M and MRS(A ) parton distribution functions. We assign an uncertainty of ±6%, 
which reflects the uncertainties above, added in quadrature. We quote as our final result for 
qq — ► tb, ib a cross section of .88 ± .05 pb at the Tevatron,f] and 10.2 ± 0.6 pb at the LHC. 

An additional source of uncertainty is the top-quark mass. A plot of the next-to-leading 
order cross section for qq — > tb, ib as a function of the top-quark mass is show in Fig. [| 
It is anticipated that this uncertainty will be ±6 GeV when the data from Run I at the 
Tevatron are fully analyzed. This yields an uncertainty of ±15% in the cross section at the 



Tevatron. The uncertainty in the mass is expected to decrease to ±4 GeV in Run II [24 



corresponding to an uncertainty in the cross section of ±10%. A high-luminosity Tevatron 
might be capable of reducing the uncertainty in the mass to ±2 GeV [23|, which would yield 
an uncertainty in the cross section of ±5%. The uncertainty in the cross section at the LHC 
is comparable. 

Much can be done to reduce the uncertainty in the calculation. The next-to-next-to- 
leading-order correction to the Drell-Yan process is already known ||. The full next-to- 
next-to-leading-order QCD correction to qq — > tb can and should be completed in the near 
future. This should reduce the uncertainty in the cross section from yet higher orders to 
below the 1% level. A reliable estimate of the uncertainty in the parton distribution functions 
requires a set with built-in uncertainties, which we hope will be available in the near future. 

It seems likely that by the time the process qq — > tb is observed in Run II at the Tevatron, 
the theoretical uncertainty in the cross section will be slightly larger than ±10%, due mostly 
to the uncertainty in the mass. This is adequate in comparison with the anticipated exper- 
imental errors. The statistical error on the measured cross section in Run II will be about 
±20% 0. This corresponds to a measurement of \ V t b\ with an accuracy of ±10% (assuming 
\Vtb\ ~ 1). A high-luminosity Tevatron, which could potentially deliver 30 fb _1 over several 

9 For y/S = 1.8 TcV, the cross section at the Tevatron is .73 ± .04 pb 
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Figure 10: Cross section for qq — > W ± H and qq — > ZH at the Tevatron and the LHC, calculated at 
next-to-leading order using CTEQ3M parton distribution functions. 



years, would allow a measurement of the cross section with a statistical uncertainty of about 
6%, with a comparable theoretical uncertainty. Combining the statistical and theoretical 
uncertainties in quadrature, this corresponds to a measurement of |V#,| with an accuracy of 
about ±4%. 

The process qq — > tb is also important as a background to the process qq —>■ WH with 
H —>■ bb at the Tevatron. We show in Fig. [10] the next-to-leading-order cross sections for 
qq — > W ± H, as well as qq — > ZH, at both the Tevatron and the LHC [25|]. The significant 
increase in the qq — > tb cross section at next-to-leading order could have a negative impact 
on the ability to find an intermediate-mass Higgs boson at the Tevatron. 
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Appendix 

Below are the form factors for the Yukawa correction to the matrix element of the tb 
charged current. These corrections arise from loops of Higgs bosons and the unphysical 
scalar W and Z bosons associated with the Higgs field in the gauges. The masses- squared 
of the unphysical scalar bosons are £M^,£Mf. In the numerical calculations, we set £ = 
(Landau gauge). The integrals were reduced to the standard one-, two-, and three-point 
scalar loop integrals and then evaluated with the aid of the code FF pfifl . The notation is 
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adopted from |TJ| ; the arguments of the functions give the internal masses-squared followed 
by the external momenta squared. 



F 1 = -[AC 2 ^M^M 2 H ,m 2 f ,q 2 ,m 2 ,0)+AC 2 ^M^^M 2 z ,m 2 ;q 2 ,m 2 ,0) 

+ B 1 (M 2 H , m 2 ; m 2 t ) + (M 2 H - 4m 2 )B' (M 2 H , m 2 ; m 2 ) 
+ B^M 2 Z , m\- m\) + ^M 2 Z B' ^M 2 Z , m 2 ; m\) 
+ B 1 (£M 2 V , 0; m\) + (£M^ - m^B^M^, 0; m 2 ) 
+ Bi{£M*y, m 2 ; 0) + {iM 2 w + m 2 t )B' (£M^, m 2 ; 0)] 



^[C 23 (eM^, M 2 H , m 2 ; q 2 , m 2 , 0) + C 23 (£M^, ^Mf, m^; q 2 , m 2 , 0) 
f C 21 (eM^, M|, mf; g 2 , m 2 , 0) + C 21 (£M^, ^Mf, m 2 ; g 2 , m 2 , 0) 
f 2Cn(eM^,M|,m 2 ;g 2 ,m 2 ,0)] 



Note: In the reduction of the three-point integrals, a misprint was discovered in Ref. [19 |. 
On p. 199 in Appendix E, in an unnumbered equation near the bottom of the page, C 22 and 
C 2 3 were transposed. The correct equation is (C 2 3,C 22 ) = X _1 (_R 4 , i? 6 ). 
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